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I.  INTRODUCTION 


Reactions  between  glass  and  an  aqueous  solution  are  important  In  a  number 
of  technical  areas.  This  is  particularly  true  for  the  proposed  encapsulation 
of  radioactive  wastes  in  vitreous  media. In  this  case  the  possible  leach¬ 
ing  by  water  of  radioactive  materials  from  the  encapsulation  matrix  is  of 
prime  concern. 3  The  liquid-glass  interaction  is  comprised  of  a  number  of 
complicated  processes;  the  importance  of  each  depends  on  glass  composition, 
reaction  conditions,  and  the  period  of  reaction.*  Studying  these  processes 
during  the  early  periods  of  glass  exposure  to  solution  results  in  some  simpli¬ 
fication.  The  "short  time  limit"  investigation  reduces  the  effects  of  surface 
etching  yet  still  provides  valuable  information  concerning  the  chemical  pro¬ 
cesses  occurring  near  the  glass  surface.^ 

Franek  et  al.^’^  have  investigated,  both  experimentally  and  theoreti¬ 
cally,  during  the  short  time  regime,  the  reactions  between  silicate  glass 
surfaces  and  aqueous  cesium  chloride  (CsCl)  solutions.  One  interesting 
process  observed  was  termed  a  "phase-boundary  reaction."  In  this  process  an 
aqueous  species,  Cs+,  is  believed  to  undergo  a  surface  reaction  prior  to  its 
diffusion  into  the  glass.  Their  experimental  data,  presented  as  penetration 
profiles  of  C820  into  the  glass  as  a  function  of  the  glass's  solution-exposure 
time,  were  analyzed  in  terms  of  a  quasi-emplrlcal  formulation.  This  model 
considered  only  a  single  species,  Cs+,  which  was  allowed  to  diffuse  into  the 
glass  subject  to  a  time-varying  surface-concentration  boundary  condition. ^ 

The  surface-concentration  boundary  condition  was  taken  empirically  from  the 
experimental  data.  This  procedure  results  in  a  closed-form  expression  for  Cs 
concentration  as  a  function  of  penetration  depth,  exposure  time,  and  Cs  diffu¬ 
sion  coefficient. 

While  fairly  good  agreement  between  the  experimental  data  and  calculated 
profiles  is  obtained,  yielding  a  diffusion  coefficient  of  5  x  10“^  cm^/s,  the 
analysis  is  unsatisfying  for  several  reasons.  In  order  to  extract  the  maximum 
chemical  and  mechanistic  information  from  the  data,  it  is  desirable  to  model 
the  entire  reaction-diffusion  process  without  empirically  introducing 


information.  The  reaction  of  aqueous  Cs*  with  glass  surface  sites  to  yield 
additional  Cs  species  that  diffuse  into  the  glass  is  inherently  a  multi- 
component  process.  Nonempirlcal  modeling  is  then  virtually  Impossible  within 
the  framework  of  a  single-component  model.  Also,  parameters  extracted  from 
the  modeling,  e.g.  the  diffusion  coefficient,  are  model-dependent.  It  would 
be  highly  desirable  to  compare  penetration  profiles  predicted  by  different 
physically  meaningful  reaction-diffusion  models  to  study  the  effects  on  the 
extracted  parameters. 

In  this  report,  the  phase-boundary  reaction  data  of  Franek  et  al.^  are 
reanalyzed  in  terms  of  two  different  mechanistic  processes.  Detailed  descrip¬ 
tions  will  be  subsequently  presented.  Briefly,  as  outlined  in  Table  1,  the 
first  process  examined  is  the  reaction  of  aqueous  Cs*  with  a  fixed  number  of 
sites  on  the  glass  surface,  yielding  a  Cs  species  that  diffuses  into  the 
glass.  The  second  process  examined  is  the  rapid  reaction  of  aqueous  Cs*  with 
a  time-dependent  number  of  surface  sites,  followed  by  Cs  diffusion  into  the 
glass.  While  the  times  considered  are  brief  enough  to  preclude  significant 
surface  removal,  they  are  sufficient  for  the  etching  process  to  generate  a 
time-dependent  number  of  reaction  sites.  In  addition  to  modeling  the  pene¬ 
tration  of  Cs  into  the  glass  by  simple  diffusion,  two  other  penetration  models 
were  considered  in  this  second  process:  first,  Cs  penetration  as  described  by 
the  Doremus  Interdiffusion  model, ®  and  second,  a  two-pathway  diffusion  mechan¬ 
ism.  The  lnterdiffuslon  and  double-diffusion  penetration  models  were  not 
developed  for  the  first  process.  The  reasons  for  this  are  given  in  the  text. 

All  of  these  models  yield  partial  differential  equations  (PDGs)  that  are 
solved  by  means  of  standard  finite-difference  techniques.  The  principal 
advantage  of  this  approach  is  that  it  allows  a  consistent  study  of  different 
physically  plausible  processes.  By  studying  the  model  solutions,  one  gains 
information  about  which  mechanisms  actually  occur  within  the  glass.  This 
paper  also  demonstrates  the  value  of  formulating  mechanistic  models  in  terms 
of  PDEs  and  then  directly  solving  these  equations.  This  technique,  which 
provides  a  means  of  treating  general  multicomponent  reaction-diffusion 
processes,  appears  not  to  have  previously  received  the  attention  its  utility 
warrants. 
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Table  1.  Phase -Boundary  Reaction  Processes 


Process 

Surface  Process 

Glass 

Penetration 

Model 

1 

2 

3 

1 

Cs+(aq)  reacts  with  fixed 
number  of  glass  surface 
sites  to  yield  diffusing 
species,  Cb(s) 

Simple 

diffusion 

— 

— 

2 

Cs+(aq)  reacts  with  time- 
dependent  number  of  glass 
surface  sites  to  yield 
diffusing  species,  Cs(s) 

Simple 

diffusion 

Inter¬ 

diffusion 

Double 

diffusion 
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II.  MATHEMATICAL  DEVELOPMENT 


In  this  section  the  two  possible  chemical  processes  for  the  phase¬ 
boundary  reaction  to  be  investigated  are  described  and  their  PDEs  are 
formulated.  The  first  process  supposes  that  aqueous  Cs,  Cs(aq),  reacts  with 
unoccupied  glass  surface  sites,  G,  to  yield  a  Cs  species,  Cs(s),  that  pene¬ 
trates  into  the  glass  via  simple  one-dimensional  diffusion.  The  initial 
concentration  of  reaction  sites  on  the  glass  surface  is  a  constant,  [G]®. 
Simple  diffusion  represents  the  most  elementary  diffusion  model.  It  would  be 
expected  to  be  valid  in  cases  where  the  penetrating  species  do  not  interchange 
with  components  of  the  glas6.  Also,  it  may  be  found  to  be  valid  in  self- 
diffusion  situations,  e.g.  where  Na  diffuses  in  a  Na  silicate  glass.  Simple 

diffusion  has  been  successfully  used  to  model  diffusion  of  various  species 
q_12 

into  glasses. 7 

Equation  (1)  shows  the  surface  reaction  considered: 


C8(aq)  ♦  G 


Cs(s) 


CD 


The  equilibrium  constant,  Ke<J,  is  the  ratio  of  the  forward  reaction-rate 
constant,  Kj,  to  the  reverse  reaction-rate  constant,  K_ j .  The  general  one- 
dimensional  diffusion  equation  for  the  concentration  of  a  species  A,  [A],  is 

3[Al  ,  3 2 [A]  , 

3t  °A  3x2 


with  Da  being  the  diffusion  coefficient  of  the  A  species. 

The  PDE  for  this  first  process  is 

.  D  lifCsCsiU  K  {Cs(aq)  1  { [G] °-  fCs(s)j]  -  K  {Cs(s)J  (3) 
"t  Zx*  1  1 
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with  D  being  the  diffusion  coefficient  of  the  penetrating  species.  A  single¬ 
boundary  condition,  the  solution  concentration  of  Cs(aq),  must  be  specified, 
while  at  time  t  ■  0,  { Cs ( s ) ]  is  equal  to  zero  for  all  positions  x  within  the 
glass.  The  solution  volume  is  assumed  to  be  large  enough  so  that  the  con¬ 
centration  of  Cs(aq)  remains  constant  during  the  time  periods  of  interest.  To 
constrain  the  reaction  in  Eq.  (1)  to  only  the  glass  surface,  and  K_^  are 
zero  everywhere  except  at  the  glass  surface,  where  x  ■  0.  This  is  justified 
because  Cs(aq)  is  present  only  at  the  surface  of  the  glass.  This  simple 
modeling  of  the  first  process  is  sufficient  to  demonstrate  its  inability  to 
explain  the  observed  experimental  data.  Consequently,  more  sophisticated 
penetration  models  for  this  process  were  not  developed. 

The  second  process  considered  supposes  that  the  number  of  reaction  sites 
on  the  glass  surface  is  time-dependent.  The  number  of  active  surface  sites 
increases  as  a  function  of  time  because  of  the  exposure  of  the  glass  to  the 
aqueous  solution.  While  the  exposure  times  of  interest  are  not  sufficient  for 
the  etching  process  to  remove  significant  layers  of  glass,  it  is  known  that 
even  brief  exposure  to  aqueous  solution  can  alter  the  glass  surface.^  Once  an 
active  site  has  been  formed,  a  Cs(aq)  rapidly  comes  in  contact  with  it,  form¬ 
ing  the  penetrating  species  Cs(s).  Consequently,  in  this  process  the  rate- 
limiting  step  in  the  formation  of  Cs(s)  is  the  generaton  of  active  sites 
rather  than,  as  in  the  previous  process,  the  reaction  of  Cs(aq)  to  form  Cs(s). 

Several  chemical  reactions  are  needed  to  describe  this  process.  The 
Increasing  surface-site  concentration  is  described  by 

K' 

G(NA)  — ->  G(A)  (4) 

with  G(NA)  representing  nonactive  surface  sites,  G(A)  being  the  reactive  glass 
sites,  and  K'  being  a  rate  constant  for  this  reaction.  The  surface  activa¬ 
tion  process  may  in  reality  be  a  complicated  series  of  reactions  depending  on 
solution  constituents.  For  simplicity,  it  is  modeled  as  a  simple  irreversible 
reaction.  The  rapid  reaction  of  Cs(aq)  is  given  by 


Cs(aq)  +  G(A)  ■■.-?>  Cs(s) 


(5) 


The  rate  of  Cs(s)  formation  at  the  surface  Is  specified  by  the  rate  equation 


d[G^-  -  KJ  {[G(NA)]°  -  [G( A) ]}  -  dl£s<s)]  (6) 

with  (G(NA)]®  being  the  Initial  number  of  nonactive  sites  or,  equivalently, 
the  maximum  possible  number  of  active  sites.  Equation  (6)  may  be  solved 
directly  to  yield 

[G(A) ] (t)  -  [G(NA)]°  [1  -  exp(-K't)]  -  [Cs(s)](x  -  0)  (7) 

with  [Cs(s)](x  *  0)  representing  the  surface  concentration  of  Cs(s). 

As  a  first  approach,  the  penetration  of  Cs(s)  into  the  glass  may  be 
modeled  by  simple  diffusion,  assuming  a  constant  diffusion  coefficient,  yield¬ 
ing  the  PDE 

3[C3^)]-  -  2 [ Cs ^-S-~  +  K'{  [G(NA)]°  -  [G(A) ]}  (8) 

3x 

with  (G(A)J  specified  by  Eq.  (7).  This  model  is  most  similar  to  the  analysis 
carried  out  be  Franek  et  al. 

Solutions  of  Eq.  (8)  were  similar  to  the  experimental  data.  Conse¬ 
quently,  more  complicated  penetration  models  were  developed  for  this  possible 
phase-boundary  reaction  process.  The  second  penetration  model  considered  for 
this  process  again  uses  Eqs.  (4)  and  (5)  for  the  surface  reaction.  However, 
the  Cs  penetration  into  the  glass  is  assumed  to  occur  by  an  interdiffusion 
mechanism.  In  this  model,  for  each  Cs  (typically  assumed  to  be  in  an  ionic 
form)  diffusing  into  the  glass,  an  ion  within  the  glass  diffuses  back  into  the 
solution.  Generally,  it  would  be  expected  that  alkali  ions  within  the  glass 
would  be  the  back-dif fusing  species. ^  The  interdiffusion  process  has 
been  particularly  effective  in  describing  the  exchange  of  ions  in  solution, 
such  as  Na+,  H*,  and  Ag+,  with  ions  in  the  glass. 5.13-16  Qne  aSpect  Qf  this 
model  is  that  the  diffusion  coefficient  depends  on  the  concentrations  of  the 
various  interdif fusing  ions  within  the  glass.  Since  these  concentrations  vary 
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with  penetration  depth,  the  diffusion  coefficient,  referred  to  as  the  inter¬ 
diffusion  coefficient,  is  position-dependent.  The  resulting  PDE  cannot  be 
solved  in  closed  form. 

For  interdlf fusion,  the  diffusion  equation  takes  the  form  (see,  for  exam¬ 
ple,  Doremus,  Ref.  8). 

(9a) 

with  D,  the  Interdiffusion  coefficient,  given  by 


(r r-  OtA]  +  W 


Da  and  Dg  are  the  self-diffusion  coefficients  of  the  penetrating  ion  A,  and 
the  back-diffusing  ion  B,  respectively.  The  total  concentration  of  A  and  B, 
which  is  equal  to  the  B  concentration  at  t  ■  0,  is  indicated  by  [T] . 

The  resulting  PDE  to  be  solved  is 


3 [ Cs ( s ) ] 

at 


a‘[cs(s)i 

„_2 


bWDcs 

(IT]  +  [Cs(s)]b)2 


a[Cs(s)I  '2 


+  R'{[G(NA))°  -  [G(A)]} 


(10a) 


»  Cs 

b.—  1 


(10b) 
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The  final  penetration  model  to  be  considered  for  this  process  treats  the 
possibility  of  glass  penetration  by  multiple  diffusion  paths.  A  number  of 
authors  have  experimentally  observed  what  appears  to  be  alkali  diffusion  via 
multiple  pathways.  Frischat1^  has  been  able  to  fit  Na  penetration  profiles  in 
several  alkali-silicate  and  alumlno-slllcate  glasses  to  models  using  two  inde¬ 
pendent  diffusion  pathways.  Similiarly,  "double-diffusion**  mechanisms  were 
believed  responsible  for  Na  penetration  profiles  observed  by  Schaeffer  and 
Mecha*®  in  high-silica-content  glasses.  Lim  and  Day*®  also  found  that  a 
multiple-diffusion  mechanism  yielded  the  best  agreement  with  their  experi¬ 
mental  Na  diffusion  profiles.  They  specifically  discussed  penetration  in 
terms  of  vacancy  and  interstitialcy  mechanisms.  In  light  of  experimental 
evidence  indicating  multiple  diffusion  pathways  for  Na  penetration,  it  was 
felt  that  8lmlliar  mechanisms  for  Cs  penetration  should  be  considered. 

To  model  a  double-diffusion  process,  two  different  active  glass  sites,  G* 
and  G^,  are  formed  via 


and 


G*(NA)  — — >  Gl(A) 


g2(na)  — ?->  G2(A) 


(11a) 


(lib) 


Each  of  the  two  types  of  active  sites  is  assumed  to  rapidly  acquire  Cs,  form¬ 
ing  two  independently  diffusing  species,  Cs*(s)  and  Cs^(s).  Two  different 
diffusion  coefficients,  Dj  and  D2,  characterize  the  diffusion  of  these 
species.  Each  of  these  diffusion  processes  has  been  treated  as  simple  dif¬ 
fusion.  For  this  penetration  model,  two  coupled  PDEs,  Eqs.  (12a)  and  (12b), 
are  obtained: 

-  D.  +  iq  { [G1(NA)J°  -  [G1(A))}  (12a) 

3t  1  3x/  1 


3-tc?-2(_8J-L  -  D,  3-2.L?8.2i8-ll  +  k;  { [G2(NA)J°  -  [G2(A) )}  (12b) 

dt  1  ax'*  1 
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It  would  have  been  straightforward  to  treat  each  as  an  Interdiffusion  mech¬ 
anism*  However,  incorporation  of  the  interdiffusion  mechanism  would  require 
the  specification  of  additional  parameters.  Considering  the  limited  quantity 
of  data  to  be  analyzed,  this  was  not  considered  practical. 

These  ?DEs  may  be  solved  approximately  by  using  standard  finite- 
difference  techniques. » 21  Briefly,  the  glass  is  treated  as  a  semi-infinite, 
one-dimensional  slab.  During  the  exposure  times  of  interest,  Cs  penetrates 
only  a  limited  distance  into  the  glass.  Consequently,  a  finite  grid  in  both 
position  x  and  exposure  time  t  may  be  used.  The  grid-point  spacing  for  posi¬ 
tion  and  time  are,  respectively,  Ax  and  At.  Within  this  grid  system  it  is 
possible  to  approximate  derivatives  by  finite  differences.  Denoting  the  value 
of  Cs  concentration  at  the  grid  point  corresponding  to  x  ■  iAx  and  t  -  jAt  by 
j,  the  simple  derivative  formulas  are 


Substitution  of  these  expressions  into  each  PDE  results  in  an  "explicit" 
system  of  algebraic  equations  for  U  at  each  grid  point.  The  system  is  ex¬ 
plicit  in  the  sense  that  if  all  values  of  j  for  a  specific  t  ■  jAt  are 
known,  values  of  each  may  be  solved  for  individually. 

While  implementation  of  these  methods  is  straightforward  for  process  2, 
certain  complications  arise  for  process  1.  In  process  1,  Cs(aq)  reacts  with  a 
fixed  number  of  surface  sites  to  yield  the  diffusing  species,  Cs(s).  In  terms 
of  the  numerical  solution,  this  process  is  assumed  to  be  composed  of  three 
discrete  substeps,  each  of  which  occurs  during  a  single  time  increment,  At. 


First,  Cs(aq)  undergoes  the  surface  reaction.  The  second  and  third  terms  on 
the  right-hand  side  of  Eq.  (3)  account  for  this,  essentially  generating  a 
time-dependent  boundary  condition  for  the  diffusion  portion  of  the  equation. 
The  second  substep  during  At  is  the  diffusion  of  the  Cs(s),  generated  at  the 
surface,  into  the  glass.  It  is  crucial  to  note  that  diffusion  of  Cs(s)  away 
from  the  surface  reduces  the  Cs(s)  concentration  at  the  surface.  (This  is  in 
contrast  to  process  2,  where  the  reaction  of  Cs(aq)  is  assumed  to  be  rapid 
enough  to  Immediately  replenish  the  surface  concentration.)  Consequently,  the 
third  substep  is  the  reduction  of  the  Cs(s)  surface  concentration  as  a  result 
of  diffusion.  The  amount  of  Cs(s)  that  has  diffused  into  the  glass  during 
At  is  the  difference  between  the  amount  prior  to  and  after  At.  This  amount 
of  Cs(s)  must  be  removed  from  the  surface  layer.  A  reasonable  approximation 
to  process  1  is  provided  by  this  scheme. 

The  power  of  these  finite-difference  techniques  lies  in  their  ability  to 
solve  directly  the  PDEs  resulting  from  complicated  mechanisms.  The  mechanisms 
analyzed  here  are  relatively  simple,  a  result  of  limited  experimental  data. 
Within  this  formalism,  though,  significantly  more  complicated  mechanisms  may 
be  studied.  It  is  a  straightforward  process  to  generalize  to  additional 
components,  interactions  between  components  within  the  glass,  more  involved 
surface  reactions,  etc. 
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III.  INTERPRETATION  OF  EXPERIMENTAL  DATA  AND  COMPUTATIONAL  PROCEDURES 


The  experiments  performed  by  Franek  et  al.^  consisted  of  immersing  glass 
samples  in  CsCl  solutions  at  90° C  for  various  periods  of  time.  Penetration 
profiles  were  obtained  by  using  a  secondary- ion,  mass-spectrometric,  neutral- 
particle-beam  technique.  To  reanalyze  the  experimental  data,  the  CS2O  concen¬ 
tration  profiles  represented  in  Fig.  1  of  their  paper  for  exposure  times  of 
10,  30,  60,  and  85  min  were  digitized  with  a  Hewlett-Packard  9825T  mini¬ 
computer.  Errors  due  to  digitization  are  less  than  the  size  of  the  symbols 
representing  the  data  points  and  much  less  than  the  experimental  errors 
Indicated  in  Fig.  2  of  their  paper. 

For  comparison  of  the  experimental  data  with  the  models  described  pre¬ 
viously  in  this  report,  the  C82O  concentration  profiles  in  Franek  et  al.'s 
paper  given  in  terms  of  molZ  were  converted  to  units  of  atoms/cm3  by  using  the 
indicated  composition  of  the  technical  glass  used.  In  addition,  the  aqueous 
Cs  concentration,  Cs(aq),  was  6  x  10^  ions/cm3,  whereas  the  total  alkali  ion 
content  [T]  in  this  glass,  a  necessary  parameter  in  the  interdif fusion  model, 
was  determined  to  be  7.2  x  lO^1  ions/cm3. 

To  solve  the  PDEs  with  the  finite-difference  technique,  solutions  were 
calculated  over  a  penetration  depth  of  25  nm,  with  a  grid  point  spacing  equal 
to  0.25  nm.  The  time  step  was  typically  1  s.  If  the  time  step  was  made  too 
large,  the  method  became  unstable  and  the  solution  rapidly  diverged.  The 
finite-difference  equations  were  solved  with  a  Hewlett-Packard  9825T  mini¬ 
computer.  Solutions  for  exposure  times  up  to  5100  s  typically  required  about 
A  h  of  computer  time. 

When  fitting  the  experimental  data  to  the  computer-generated  solutions, 
there  are  several  free  parameters  to  be  specified.  For  process  1,  values  are 
needed  for  the  surface  reaction  rate  constants  Kj  and  ,  the  surface  site 
concentration  [G]^,  and  the  diffusion  coefficient  D.  The  surface- reaction 
rate  constants  specify  how  rapidly  Cs  builds  up  on  the  glass  surface,  whereas 
the  diffusion  constant  describes  how  rapidly  the  Cs  penetrates  the  glass.  The 
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fits  were  found  to  be  sensitive  to  values  of  Kj  but  relatively  insensitive  to 
This  is  due  to  the  short  tine  regimes  studied.  Over  these  time  periods, 
insufficient  Cs(s)  Is  formed  to  make  the  reverse  reaction  rate  very  signifi¬ 
cant.  Consequently,  the  fits  were  not  very  sensitive  to  variations  in  the 
ratio  of  Kj  to  the  surface-reaction  equilibrium  constant,  when  K_j  was 

changed.  Typically,  Keq  was  set  at  5  x  lO"1^  cm^/atom,  indicating  that  at 
complete  equilibrium  most  surface  sites  would  be  filled.  Ke(j  could  be  varied 
several  orders  of  magnitude  by  changing  K_j  without  disturbing  the  fit.  [G]® 
and  Kj  are  found  to  be  the  important  parameters  in  attempts  to  fit  the  surface 
concentration.  Adjustment  of  D  varies  the  penetration  profile. 

In  the  case  of  process  2  with  a  simple-diffusion  penetration  model,  the 
site-activation  rate  constant  ft',  the  maximum  surface-site  concentration, 
(G(NA)]0,  and  D  must  be  specified.  For  this  process,  the  important  parameters 
for  mdellng  the  surface  concentration  are  found  to  be  [G(NA)]®  an<j  By 

varying  (G(NA)]^,  one  produces  a  range  of  values  resulting  in  good  fits, 
although  for  a  given  [G(NA)]®  only  a  single  ft'  yields  the  best  fit.  This 
ambiguity  might  have  been  resolved  if  additional  data  had  been  taken  for 
somewhat  longer  periods  of  time.  As  exposure  time  increases,  the  Cs  concen¬ 
tration  at  the  glass  surface  begins  to  saturate.  Observation  of  this  satur¬ 
ation  would  allow  unambiguous  determination  of  K'  and  [G(NA) ]®.  Our  criterion 
was  to  find  the  minimum  surface-site  concentration  for  which  could  be 
adjusted  to  yield  a  good  fit.  Although  this  ambiguity  exists,  it  does  not 
affect  our  results  with  regard  to  the  most  significant  surface-boundary 
reaction  model.  The  diffusion  coefficients  were  easily  adjusted  to  match  the 
experimental  penetration  profiles.  The  parameter  set  obtained  for  the  simple- 
diffusion  penetration  model  formed  the  core  of  the  parameter  sets  used  in 
subsequent  penetration  modeling.  Generally,  ten  computer  solutions  were 
required  to  optimize  the  parameter  set  for  each  model. 
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IV.  RESULTS  AND  DISCUSSION 


Calculations  based  on  process  1  were  unable  to  match  the  experimental 
results  while  maintaining  physically  meaningful  parameter  values.  To  obtain 
reasonable  agreement  between  the  experimental  data  and  calculation  results 
required  a  maximum  active-site  concentration  of  [G(NA)]®  >  1  x  10^  sites/cm^. 
Since  a  typical  glass  atom  density  would  be  approximately  1  x  10^  atoms/cm'*, 
the  required  density  was  considered  unreasonable.  Additionally,  even  when 
unrealistic  parameter  sets  are  employed,  process  1  results  do  not  agree  with 
the  data  as  well  as  process  2  results  using  realistic  parameters.  These 
calculations  indicate  that  it  is  highly  unlikely  that  the  phase-boundary 
reaction  proceeds  via  a  klnetically  limited  reaction  of  Cs(aq)  with  active 
surface  sites  on  the  glass. 

In  contrast,  results  from  process  2,  the  rapid  reaction  of  Cs(aq)  with  a 
time-dependent  number  of  surface  sites,  match  experimental  data  quite  well. 

The  general  quality  of  the  fits  allows  detailed  modeling  of  the  penetration 
process.  Figures  1  through  3  display  the  best  fits  to  Franek  et  al.’s  experi¬ 
mental  data  resulting  from  each  proposed  penetration  model.  Results  for 
reaction  followed  by  simple  diffusion  are  presented  in  Fig.  1.  The  diffusion 
coefficient  was  found  to  be  4.5  x  10“*7  cm^/s.  The  current  result  is  essen¬ 
tially  the  same  as  Franek  et  al.’s7  value  for  the  diffusion  coefficient, 

5  x  10" 17  cm^/s.  This  is  reasonable,  as  the  current  results  are  essentially 
the  exact  solution  of  the  model  previously  treated  empirically. 

Reaction  followed  by  an  lnterdlffusion  process  is  expected  to  be  a  more 
physically  realistic  model.  Two  regimes  of  the  lnterdlffusion  process  were 
investigated,  Dg  >  Dgg  and  Dg  <  Dgs.  The  best  lnterdlffusion  fits  were  found 
in  the  range  5  Dgg  >  Dg  >  0.2  Dgg,  with  Dg  -  Dqs  yielding  the  simple  diffusion 
fit.  In  Fig.  2  we  present  representative  fits  for  both  lnterdlffusion 
regimes.  For  Dg  -  0.2  DCg,  DCg  was  found  to  be  1  x  10”16  cm^/s,  whereas  for 
Dg-  5  D^g,  DCg  was  found  to  be  3.4  x  10”*7  cmfys.  Each  regime  represents  a 
different  physical  situation.  The  significance  of  these  regimes  will  be 
subsequently  discussed. 
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rlg.  1.  Cs  Penetration  Profiles  from  the  Slnple  Diffusion 

Penetration  Model  (solid  lines)  and  Franek  et  al.'s 
Experimental  Data  for  Exposure  Times  of  600  (•), 
1800  (■),  3600  (O),  and  5100  (□)  s 


[Cs]  latoms/cnr)  [Cs]  (atoms/crrrl 


Fig.  2.  Cs  Penetration  Profiles  from  Interdiffusion 
Penetration  Models,  (a)  Dg  -  5  DCg  and 
(b)  D®  •  0.2  DCg  are  given  by  solid  lines. 
Experimental  data  are  as  in  Fig.  1. 
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Results  for  the  third  penetration  model,  reaction  followed  by  double 
diffusion,  are  shown  In  Fig.  3.  The  penetration  profiles  of  Cs  diffusing  by 
each  pathway  are  shown  along  with  the  total  Cs  penetration  profile.  The 
experimental  data  were  best  described  as  the  sum  of  the  major  process 
(Dj  ■  3  x  10" 17  cm^/s)  and  the  minor  process  (D2  ■  2  x  10“^  cm^/s).  It  Is 
interesting  to  note  that  to  obtain  the  best  fit  to  the  data,  the  minor  pene¬ 
tration  pathway  was  required  to  have  Its  surface  concentration  saturated  at  an 
exposure  time  of  ~3600  s.  Consequently,  a  unique  pair  of  values  could  be 
found  for  K'  and  [G^(NA)]®,  4.2  x  10“^  cm^/atom  •  s  and  4  x  10^®  atoms /cm^, 
respectively.  Unfortunately,  as  noted  previously,  the  main  pathway  was 
unsaturated  for  all  exposure  times  considered,  which  prevented  a  unique  deter¬ 
mination  of  KJ  and  [G*(NA)]®. 

Can  one  penetration  model  be  said  to  describe  the  experimental  results 
most  accurately?  Each  model  yields  reasonable  agreement  with  the  data  for 
penetration  depths  less  than  5  nm.  However,  beyond  5-nm  penetration  depths 
the  fits  differ  significantly.  In  Fig.  4  the  penetration  profiles  beyond  5  nm 
for  each  model  are  displayed  for  comparison.  Inspection  of  Fig.  4  Indicates 
that  the  double-diffusion  mechanism  yields  the  best  fit  to  the  experimental 
data.  To  quantify  the  differences  in  quality  of  fits,  an  average  absolute 
percent  difference.  A,  between  the  experimental  point  and  the  predicted 
value  has  been  computed  for  each  model.  After  consideration  of  all  the 
experimental  data,  A’s  for  double  diffusion,  interdiffusion  (Dg  >  Dqs),  simple 
diffusion,  and  lnterdlffusion  (Dg  <  Dq8)  are  respectively  20%,  27%,  50%,  and 
59%.  If  only  data  beyond  a  5-nm  penetration  is  considered,  the  respective 
A's  are  25%,  39%,  82%,  and  96%.  The  double-diffusion  mechanism  shows  the  best 
agreement  with  both  the  entire  set  of  experimental  results  and,  In  particular, 
with  the  experimental  data  beyond  5  nm  in  depth.  Consequently,  of  the  models 
currently  studied,  the  double-diffusion  model  appears  to  explain  most 
accurately  the  penetration  processes  occurring  within  the  glass.  Of  course, 
this  conclusion  might  not  be  valid  should  any  significant  errors  have  occurred 
during  the  acquisition  of  the  experimental  data. 


Cs  Penetration  Profiles  from  (a)  Simple  Diffusion, 
(b)  Interdiffusion  (D»  *  5  Dgs),  (c)  Interdiffusion 
(Db  *  0.2  Dgs),  and  (a)  Double-Diffusion  Penetration 
Models  over  Penetration  Depths  of  5  to  25  nm. 
Experimental  data  are  as  in  Fig.  1. 


A  side  aspect  is  that  the  model  providing  the  second-best  fit  is  inter- 
diffusion,  with  Djj  *  0.2  Dqs.  The  most  reasonable  Interdiffusion  process  is 
the  replacement  of  alkali  ions  in  the  glass  with  Cs  ions  from  the  solution. 

The  alkali  ions  in  the  glass  are  Na+  and  1C*-,  which  would  be  expected  to 
diffuse  more  rapidly  than  Cs+.^  Consequently,  this  model,  even  though  it 
yields  reasonable  agreement  with  experiment,  should  most  likely  be  discounted. 
A  possible  reason  for  the  inability  of  the  interdiffusion  mechanism  to  realis¬ 
tically  describe  the  Cs  penetration  process  is  that  Dqs  and  Dg  may  be 
concentration-dependent.  Varshneya  and  Milberg^*^  have  shown  that  the 
potassium  self-diffusion  coefficient  for  temperatures  >  340°C  can  be 
concentration-dependent.  If  experimental  self-diffusion  coefficient  data  for 
the  ions,  glass  type,  and  temperature  under  study  existed,  it  would  have  been 
straightforward  to  incorporate  them  expicitly  into  the  finite-difference 
model.  In  the  absence  of  these  data,  these  coefficients  were  treated  as 
concentration-independent . 

With  the  observation  of  a  double-diffusion  process,  it  would  be  desirable 
to  identify  what  microscopic  mechanisms  are  responsible  for  each  diffusion 
pathway.  The  present  studies  do  not  directly  address  this  point.  However,  in 
the  light  of  previous  studies,  some  speculation  is  possible.  Lim  and  Day^ 
found  that  two  principal  diffusion  mechanisms  occur  when  Na  diffuses  into  Na, 
K,  or  Rb  silicate  glasses.  The  first  is  the  vacancy  mechanism,  in  which  a 
single  penetrating  ion  jumps  from  lattice  vacancy  to  lattice  vacancy. 
Interstitialcy,  the  second  mechanism,  is  a  two-atom  process.  A  penetrating 
ion  initially  in  an  interstitial  site  jumps  to  a  lattice  site,  displacing  the 
ion  at  that  site  to  another  interstitial  site.  In  their  experiments,  Lim  and 
Day  found  that  the  interstitialcy  mechanism  dominates  at  low  temperatures  (T  < 
300°C),  whereas  at  higher  temperatures  the  ions  diffused  via  the  vacancy 
mechanism.  Although  the  effects  of  glass  composition  and  differing  experi¬ 
mental  conditions  on  the  diffusion  mechanisms  are  not  clear,  a  simple  extrapo¬ 
lation  of  these  results  to  Franek  et  al.'s  conditions  would  indicate  that  the 
major  diffusion  mechanism  would  be  interstitialcy.  Consequently,  a  reasonable 
speculation  is  that  this  mechanism  is  the  major  diffusion  mechanism  observed. 
The  minor  mechanism  would  similiarly  be  expected  to  be  the  vacancy  mechanism. 


Finally,  we  note  that  to  really  test  a  particular  phase-boundary  reaction 
and  penetration  mechanism,  additional  experimental  data  are  required.  As  was 
stated  previously,  somewhat  longer  exposure  times  would  allow  unique  deter¬ 
mination  of  the  surface-site  rate  constant  and  the  glass-surface-site  concen¬ 
tration  used  in  the  present  modeling.  Additionally,  performing  experiments  at 
a  variety  of  aqueous  Cs  concentrations  would  allow  more  detailed  investigation 
of  the  surface  reactions.  More  complicated  reaction  mechanisms  could  then  be 
formulated  and  tested. 
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V.  CONCLUSIONS 


This  report  demonstrates  how  effectively  model  processes  may  be  compared 
to  experimental  data  through  the  direct  solution  of  the  resulting  PDEs  by 
using  a  simple  finite-difference  technique.  The  capability  to  test  a  variety 
of  model  processes  is  Important  for  several  reasons.  By  investigating  differ¬ 
ent  models,  one  gains  Insights  into  the  actual  processes  occurring  within  the 
glass,  Including  microscopic  aspects  of  the  reaction  and  subsequent  penetra¬ 
tion.  Additionally,  parameters  extracted  from  the  experimental  data,  e.g. 
diffusion  coefficients,  were  found  to  be  model-dependent;  consequently,  to 
obtain  their  proper  values,  one  must  use  the  correct  model.  Specifically, 
processes  differing  in  surface  reactions  have  been  applied  to  the  phase¬ 
boundary  reaction  data  of  Franek  et  al.  The  experimental  data  are  entirely 
consistent  with  the  rapid  reaction  of  aqueous  Cs+  with  a  time-dependent  number 
of  surface  sites.  The  required  increase  in  surface-site  concentration  with 
time  is  believed  to  occur  through  reactions  between  the  glass  surface  and  the 
aqueous  solution.  Subsequent  penetration  of  Cs  Into  the  glass  was  more 
accurately  modeled  by  a  double-diffusion  mechanism. 
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LABORATORY  OPEEATIOHS 


The  laboratory  Operation*  of  The  Aerospace  Corporation  la  conducting 
experimental  and  theoretical  Investigation  neceasary  for  the  evaluation  and 
application  of  scientific  advsnees  to  new  military  apace  systems.  Veraetlllty 
and  flexibility  have  been  developed  to  a  high  degree  by  the  laboretory  person¬ 
nel  In  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  systema.  Expertise  In  the  latest  scientific  developments  Is 
vitsl  to  the  accomplishment  of  tasks  related  to  these  problems.  The  labora¬ 
tories  that  contribute  to  this  resesrch  are: 

Aerophyslcs  laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  environmental  hasarde,  trace  detection;  spacecraft  structural 
mechanics,  contamination,  thermal  and  structural  control;  high  temperature 
thermomechanlcs,  gas  kinetics  and  radiation;  cw  and  pulsed  laser  development 
Including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control, 
atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optica,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emleslon,  photosensitive  materials  and  detectors,  atomic  frequency 
standards,  and  environmental  chemistry. 

Computer  Science  laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  space bo me 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence  and 
microelectronics  applications. 

Electronics  Research  laboratory:  Microelectronics,  GaAs  low  noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenoawna,  quantum  electronics,  laser  communications,  Ildar,  and  electro¬ 
optics;  communication  sciences,  applied  electronics,  semiconductor  crystal  and 
device  physics,  radiometric  Imaging;  millimeter  wave,  microwave  technology, 
and  RF  systems  research. 

Materials  Sciences  laboratory:  Development  of  new  isaterlals:  metal 
matrix  composites,  polymers,  and  new  forms  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  reliability;  fracture  mechanics  and 
stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures  as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  laboratory:  Magnetos pherlc,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  Interactions,  aagnetospherlc  plasma  waves;  atmospheric  and 
Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote 
sensing  using  atmospheric  radiation;  solar  physics.  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  redlatlons  on  space  system;  specs 
Instrumentation. 


